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We have performed CDCC calculations for the 6 Li + 59 Co, 144 Sm and 208 Pb systems, to investigate 
the dependence of the relative importance of nuclear and Coulomb breakup on the target charge 
(mass) at near barrier energies. The calculations were in good agreement with the experimental 
elastic scattering angular distributions for these systems and then, their predictions to the nuclear, 
Coulomb and total breakup were investigated. Although the relative importance of the nuclear 
breakup is, as expected, larger for lighter targets, this effect is not very pronounced. We also 
investigate a scaling of the nuclear breakup with the target mass and we compare the predictions 
for the integrated total breakup cross sections with experimental fusion cross sections at similar 
energies. 
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I. INTRODUCTION 



The breakup of weakly bound nuclei has been a sub- 
ject of great interest in the last years PQ, both theoret- 
ically and experimentally. Among the main investiga- 
tions within this field, there are studies of the breakup 
cross section and the influence of the breakup process on 
other reaction channels. The conclusions from those in- 
vestigations may change with the characteristics of the 
weakly bound nuclei (from now on we will assume that 
they are the projectiles), the targets and energy regime 
involved. For example, halo nuclei which have extremely 
low breakup energy threshold may behave differently 
from stable weakly bound nuclei like 6 Li, 7 Li and 9 Be. 
Usually, it is assumed that Coulomb breakup predom- 
inates over nuclear breakup when heavy targets are in- 
volved. The situation may be different in the case of light 
targets. Thus, the nature of the breakup process should 
depend strongly on the target mass. 

Regarding the influence of the breakup process on the 
fusion cross sections, at the present there is a general 
qualitative understanding that breakup enhances fusion 
at sub-barrier energies, whereas it produces some sup- 
pression above the barrier [31 [3] . Concerning the energy 
dependence of the optical potential in the scattering of 
weakly bound systems, several works show a behavior dif- 
ferent from the one found for tightly bound systems. This 
behaviour is usually called the breakup threshold anomaly 
(BTA) [4]. Recently it has been shown experimentally 
that transfer processes followed by breakup may predom- 
inate over direct breakup of stable weakly bound nuclei 



at sub-barrier energies [SHE]. 

In the present work we investigate the breakup process 
evaluating separate contributions from the Coulomb and 
from the nuclear fields, as well as the Coulomb-nuclear 
interference. We perform calculations for collisions of 6 Li 
projectiles with 59 Co, 144 Sm and 208 Pb, at near-barrier 
energies. Since the target's atomic numbers are 27, 62 
and 82, one may consider 6 Li + 59 Co, 6 Li + 144 Sm and 
6 Li + 208 Pb as a light system, a medium system and a 
heavy system, respectively. The choice of these systems 
was based on the availability of elastic scattering data 
in the literature. In this way, we could check the relia- 
bility of our calculations through comparisons with the 
scattering data. 

The present paper is organised as follows. In section 
II we present the theoretical method used in the calcula- 
tions. In sections III to V we show the results and discuss 
them. Finally, in section VI we present the conclusions 
of our work. 



II. CDCC CALCULATIONS 

It is widely accepted that the most suitable approach 
to deal with the breakup process, which feeds states in 
the continuum, is the so called continuum discretized cou- 
pled channel (CDCC) method [HUTU]. In this type of 
calculations, the continuum wave functions are grouped 
in bins or wave packets that can be treated similarly to 
the usual bound inelastic states, since they are described 
by square-integrable wave functions. In the present work 
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we use the same assumptions and methodology adopted 
in the CDCC calculations of Ref . [11] , where the elastic 
scattering of the 6 Li 4- 144 Sm system was investigated. 
We present below a summary of the main points of the 
CDCC method. Further details of the calculations can 
be found in Refs. [9HII]. 

Collisions involving 6 Li projectiles are influenced by 
the continuum states, representing mainly the breakup 
of 6 Li into a deuteron and an alpha particle. Owing to 
the low threshold of this breakup reaction (S a = 1.47 
MeV), the breakup channel is strongly coupled with the 
elastic channel. Thus, it is necessary to include the con- 
tinuum in the coupled channel calculation and this can 
be done with the CDCC method. For this purpose, we 
use the cluster model in which 6 Li is described as a bound 
state of the d + a system and the breakup channel is rep- 
resented by the continuum states of this system. This 
model has been successfully used in previous CDCC cal- 
culations in collisions of 6 Li projectiles [T31 H3] • The nu- 
merical calculations were performed using the computer 
code FRESCO [14 . In the cluster model, the projectile- 
nucleus interaction is written as 



tive 3~ (E* = 2.6145 MeV) and 5" (E* = 3.1977 MeV) 
states. The deformation parameters were taken from Ref. 
[2D] . As for the 59 Co target the quintuplet of identified 
states associated with the 2 + collective excitations were 
approximated by a single level with the energy equal to 
the centroid of the multiplet and the deformation length 
corresponding to that of the combined states [3T] . 



III. RESULTS OF THE CDCC CALCULATIONS 
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V(R, r, = V Q _ T (R, r, + V d _ T (R, r, £) (1) 

where R is the vector joining the projectile's and target's 
centers, r is the relative vector between the two clusters 
(d and a), and Q stands for any other intrinsic coordinate 
describing the projectile-target system. 

In our calculations, the continuum states of 6 Li are 
discretised as in Refs. [TTJ [TS], [TB] . Thus, we do not repeat 
the details here. The interaction between the d and the a 
clusters within 6 Li is given by a Woods-Saxon potential, 
with the same parameters as in Refs. [TTJ [T51 ITD] . 

The real parts of the V Q _ T (R, r,£) and Vd_ T (R, r,£) 
interactions were given by the double-folding Sao Paulo 
potential [17]. We have assumed that the mass densi- 
ties of the d and a clusters, required for the double- 
folding calculation, can be approximated by the charge 
densities multiplied by two, whereas the mass densities 
of the 59 Co, 144 Sm and 208 Pb targets were taken from 
the systematic study of Ref. |17j . The imaginary parts of 
V a _ T (R, r,^) and Vrf_ T (R, r,£) were chosen as to repre- 
sent short-range fusion absorption. They were given by 
Woods-Saxon functions with depth W — —50 MeV, ra- 
dial parameter = 1.06 fm and diffusivity = 0.2 fm. 
These imaginary potentials correspond to taking ingoing 
wave boundary conditions. 

The present CDCC calculations include also inelas- 
tic channels, corresponding to collective excitations of 
the targets. These channels were selected according to 
the specific nuclear structure properties of the target. 
For 144 Sm, the excitations included were the one-phonon 
quadrupole (2+, E* = 1.660 MeV) and octupole (3+, 
E* = 1.8102 MeV) first order vibrations. The values of 
the deformation parameters were obtained from Ref. [18] 
and [T!5] for the quadrupole and octupole deformations, 
respectively. For the 208 Pb target, we consider collec- 
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FIG. 1: Elastic scattering angular distributions for the three 
systems investigated. The data for the 59 Co (panel (a)), 
144 Sm (panel (b)) and 208 Pb (panel (c)) targets are respec- 
tively from Refs. [22], [23] and [21]. 



In this section, we present the results of our CDCC cal- 
culations for the breakup of 6 Li projectiles, in collisions 
with 59 Co, 144 Sm and 208 Pb targets. As a preliminary 
step, we checked that the predictions of the method for 
the elastic scattering of these systems were in good agree- 
ment with the available data. The situation is illustrated 
in figure [T] which compares theory and experiment at one 
near-barrier energy for each system. The good agreement 
obtained in the case of elastic scattering justifies the use 
of the model in calculations of other processes. 
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FIG. 2: (Color online) Breakup angular distributions at 3 
near-barrier energies for the systems discussed in the text. In 
each case, we show the total breakup cross section (solid lines) 
together with the cross sections for pure Coulomb (short- 
dashed lines) and pure nuclear (long-dashed lines) breakup. 
As in the previous figure, panels (a), (b) and (c) correspond 
respectively to results for the 59 Co, 144 Sm and 208 Pb targets. 



A. Angular distributions in Li breakup 

We now turn to the calculations of breakup cross sec- 
tions, with the purpose of assessing the relative impor- 
tance of the Coulomb and the nuclear contributions to 
the breakup process. We begin with a study of angu- 
lar distributions. In figure [2] we show the breakup cross 
sections at 3 near-barrier energies for each of the sys- 
tems mentioned above. The figures exhibit the total 
breakup cross section, together with the separate con- 
tributions from Coulomb and nuclear breakup. The sys- 
tems and the collision energies are given inside each sub- 
figure. As a reference, we mention that the heights of 
the Coulomb barrier calculated with the Sao Paulo po- 
tential are V B ( 59 Co) = 12.8 MeV, V B ( 144 Sm) = 23.8 
MeV and V B ( 208 Pb) = 29.4 MeV. Note that for each 
system the lowest energy in the figure is below the bar- 



rier whereas the highest is above it. Inspecting figure 
[2] we note that at the lowest energies the Coulomb con- 
tribution tends to dominate. In fact, for energies much 
below the barrier one expects that only Coulomb breakup 
can take place, owing to the short-range of the nuclear 
forces. On the other hand, at energies above the bar- 
rier Coulomb breakup dominates in the breakup at for- 
ward angles whereas nuclear breakup tends to dominate 
at large angles. This is not surprising, since large angle 
scattering corresponds to small impact parameters, for 
which the trajectories reach small projectile-target sepa- 
rations, within the reach of the nuclear forces. 

The transition from the Coulomb dominated angular 
region to the nuclear dominated one occurs at some cross- 
ing angle, 9q, which is a function of the collision energy. 
We have investigated the energy dependence of this angle 
for the three systems of figure [2j The results are given 
in figure [3j which shows the crossing angle as a function 
of the collision energy, normalised with respect to the 
barrier height. First, we notice that the crossing angle 
grows monotonically as the collision energy decreases. In 
this way, there is a critical energy below which Coulomb 
breakup dominates at all angles. Thus, in this energy 
region there is no crossing angle. The second interesting 
feature of this figure is that the results have a very weak 
dependence on the target. Thus, 8 can be approximated 
as a single function of E/V B , even for targets in different 
mass ranges. 
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FIG. 3: (Color online) The crossing angle, 9q, above which 
nuclear breakup dominates the breakup angular distribution. 
For details see the text. 



B. Integrated breakup cross sections 

We now investigate the total breakup cross sections, 
integrating the angular distributions of the previous sub- 
section over angles. Results for the systems of the previ- 
ous figures at three collision energies are given in table [I] 
The three energies for the different systems correspond 
to approximately the same values of E c m /V B . 

Table |T] shows some interesting properties of the 
breakup cross sections. First, one notices that, as ex- 
pected, the Coulomb breakup cross sections are system- 
atically larger for heavier systems. On the other hand, 
the strongest nuclear breakup occurs for the intermediate 
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TABLE I: Integrated breakup cross section for the systems 
discussed in the text, for three collision energies. The energies 
are given in MeV and the cross sections in mb. 



A more systematic study of the relative importance 
of the Coulomb and the nuclear forces in the breakup 
process is presented in figure |4j which shows the ra- 
tio cr^" u /cr^" c , as a function of E/V B . One can observe 
that near the Coulomb barrier this ratio is larger than 
unity for the three systems. However, the relative im- 
portance of the Coulomb breakup decreases as the energy 
increases. In the case of the lightest target, the ratio be- 
comes smaller than one for E/V B > 1.1. As expected, for 
the same values of E/V B this ratio increases with the tar- 
get's charge. However, the difference between the Co and 
Sm targets, at the same relative energies, is very small. 

IV. SCALING OF THE BREAKUP CROSS 
SECTIONS 

An interesting feature of the nuclear breakup cross 
section is the so-called "scaling" law, which says that 
this quantity at high enough energies (several tens of 
MeV/nucleon) should depend on the mass number of the 
target nucleus as [25] [26] . 



mass 1 Sm target. This is probably associated with the 
distance of closest approach and the width of the effective 
barrier in this collision but a quantitative interpretation 
of this behaviour would require further study. An impor- 
tant remark is that adding the nuclear with the Coulomb 
breakup cross sections, one does not get the total breakup 
cross section. For example, in the case of 6 Li breakup 
with the 208 Pb target at 29 MeV, the Coulomb breakup 
cross section alone is larger than the total breakup cross 
section. The same happens for the other systems at the 
lowest collision energy (below the barrier). This is a con- 
sequence of destructive Coulomb-nuclear interference in 
the breakup process. Note that this effect can also be ob- 
served in the angular distributions (e.g. figure [5] for the 
208 Pb target at 29 MeV, around 6 ~ 120 degrees). The 
last column of table [i] shows the ratio (cr t B " - crj^,) / a c^- 
Note that this ratio is always less than one. This fact 
clearly shows the destructive character of the nuclear- 
Coulomb interference. 
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FIG. 4: (Color online) Coulomb to nuclear ratio of integrated 
breakup cross section for the three systems under investiga- 
tion. 



<Cc = P l +P 2 Al /3 , (2) 

where the parameters P t and P 2 depend on the projectile, 
the structure of the target and the bombarding energy. 
Ref. [35] gives results of detailed CDCC calculations of 
elastic scattering and breakup cross sections for halo 
and non-halo systems. They study collisions of 8 B (one 
proton-halo nucleus), 11 Be (one neutron-halo nucleus) 
and 7 Be (normal, non-halo nucleus) projectiles, from 
several target nuclei. It was found that the above scaling 
law works best for the nuclear breakup of the normal 
7 Be projectile. For the halo projectiles, the scaling law 
works only approximately. In fact, the cross sections for 
these projectiles show a maximum for targets of inter- 
mediate mass, decreasing for heavy targets such as 208 Pb. 

In the present work we test the scaling law at lower en- 
ergies, close to the Coulomb barrier. Plotting the nuclear 
breakup cross section as a function of A 1 / 3 , for the same 
bombarding energy, as it was done for high energies, one 
finds that they decreases with the target mass. This is 
because the cross section changes rapidly as the collision 
energy reaches the Coulomb barrier. To eliminate this 
effect, we normalise the collision energy with respect to 
the Coulomb barrier. That is, we compare cross sections 
for the same value of E cm /V B . 

In figure 5, we show the nuclear breakup cross sections 
of 6 Li as functions of A 1 / 3 . The results are for E c , m jV B — 
0.84 (panel (a)), 1.00 (panel (b)), and 1.07 (panel (c)). 
The general behavior resembles the high energy results 
of Ref. [55], for the non-halo weakly bound 7 Be nucleus. 
In fact the almost straight lines that represent the 
curves for E c . m ./V B = 0.84, 1.00, and 1.07 are fitted 
with P 1 = -14.76, -62.60 and -49.89 mb, and P 2 = 4.11, 
16.94 and 15.41 mb, respectively. The rather large and 
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FIG. 5: (Color online) Scaling of the nuclear breakup cross 

1/3 

sections as a function of A T for the systems discussed in the 
text, at Ecuu/Va = 0.84 (panel (a)), 1.00 (panel (b)) and 1.07 
(panel (c)). 




FIG. 6: (Color online) Scaling of the Coulomb breakup cross 
sections as a function of Z T for the systems discussed in the 
text. Panels (a), (b) and (c) correspond respectively to results 
at E c . m ./V B = 0.84 1.00 and 1.07. 



negative values of P 1 are presumably traced back to 
barrier penetration effects, which limit the use of the 
geometrical picture behind the scaling law. On the other 
hand, at above-barrier energies the values of P 2 , which 
supply the slopes of the curves, are practically equal. 
Accordingly, the modified scaling law presented here 
should supply a useful and easy way to estimate the 
nuclear breakup cross section at other energies close to 
the barrier, and for other target nuclei. 

One can also derive a scaling law for the Coulomb 
breakup cross section. In figure 6, we plot <r^" u vs. 
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FIG. 7: (Color online) Comparisons of the fusion and the 
breakup cross sections for collisions of 6 Li projectiles with 
the 59 Co (panel (a)), 144 Sm (panel (b)) and 208 Pb (panel (c)) 
targets, at near-barrier energies. 



Z T , for the three systems discussed in the text. Pan- 
els (a), (b) and (c) correspond respectively to results for 
E/V B = 0.84, 1.00 and 1.07. The figures show that the 
Coulomb breakup cross sections depend linearly on Z T , 
to a very good approximation. This behaviour can be 
qualitatively explained as follows. First, we point out 
that the electromagnetic coupling matrix-elements are 
proportional to Z T , which leads to a dependence in 
the Coulomb breakup cross section. On the other hand, 
the cross sections for reaction channels are proportional 
to a l/E CAn , factor [27]. Since in each panel the collision 
energy corresponds to the same E c _ m jV B ratio, and V B is 
roughly proportional to Z T , one gets a 1/Z T factor. The 
combination of the two arguments presented above leads 
to the linear dependence obtained in figure 6. 



V. COMPARISON BETWEEN FUSION AND 
BREAKUP CROSS SECTIONS 

Now we compare the predicted integrated breakup 
cross sections with the experimental values of the cor- 
responding fusion cross section. We consider the same 
systems of the previous section, and use the fusion data 
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ofRefs. 28-30^. Figure 7 shows these comparisons. First, 
one notices that the breakup cross section for the light 6 Li 
+ 59 Co system is nearly two orders of magnitude smaller 
than the fusion cross section. Another interesting point 
is that the ratio of these cross sections is roughly con- 
stant in the whole energy interval of the figure. For the 
medium-mass and heavy systems, 6 Li + 144 Sm and 6 Li 
_l_ 208pj^ ^he sit ua ti on i s different. The breakup cross 
section is dominant at low energies whereas the fusion 
cross section dominates at high energies. The transition 
takes place in the neighbourhood of the Coulomb bar- 
rier. This behaviour is due to the Coulomb contribution 
to breakup. As the energy decreases below the barrier, 
the fusion process falls off exponentially. This is because 
fusion takes place when the projectile tunnels through 
the potential barrier. On the other hand, the decrease of 
Coulomb breakup is much slower. As the collision energy 
decreases, the classical turning point increases. Thus, the 
Coulomb part of the breakup coupling becomes weaker. 
However, the electromagnetic interactions at long dis- 
tances fall as r~( A+1 ) (where A is the multipolarity of the 
interaction), which is much slower than the exponential 
fall off. 



VI. CONCLUSIONS 

We have investigated the nature of 6 Li breakup in near- 
barrier energy collisions, with targets in different mass 
ranges ( 50 Co, 144 Sm and 208 Pb). Our theoretical cross 
sections are based on CDCC calculations, which were 
shown to lead to accurate predictions of the available 
scattering data. 

For each system, we studied the importance of contri- 
butions from the nuclear and from the Coulomb fields. 



We found that at low enough energies (E c . m . < 0.9 V B ), 
the breakup process is mainly due to Coulomb forces, 
for any scattering angle. In this energy region, the nu- 
clear and the Coulomb amplitudes for the breakup pro- 
cess interfere destructively. In this way, the cross sec- 
tion for pure Coulomb breakup may be larger than the 
cross section arising from the simultaneous action of the 
Coulomb and the nuclear fields. Above the Coulomb bar- 
rier, Coulomb breakup was shown to dominate at forward 
angles whereas nuclear breakup is dominant at larger an- 
gles. The transition takes place at some angle 6* , which 
increases with decreasing energies. A study of the energy 
dependence of the transition angle indicated that it be- 
comes nearly system independent if plotted as a function 
of the energy normalised with respect to the Coulomb 
barrier. 

Wc have shown that the nuclear breakup cross section 
has a nearly linear dependence of A 1 / 3 as suggested in 
Ref. for collisions at higher energies. On the other 
hand, the Coulomb breakup cross section was shown to 
depend linearly of the target charge, as could be pre- 
dicted by qualitative arguments. 

Finally, we made a comparison of our theoretical 
breakup cross sections with the available fusion data. 
We concluded that the latter are about two orders of 
magnitude larger than the former for the light 6 Li + 
50 Co system. For heavier systems, the breakup cross sec- 
tion becomes more important below the Coulomb barrier. 
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